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IV. C. 3  Errors  in  jump  in  velocity  versos  N/S 


ALGORITHM 


OF  THE 

MODIFIED  PROGRAM  GE02D 


If  IENTER  -  1,  then  : 

creation  of  time  traces/ 
adding  noise  if  desired/ 

If  IENTER  <  2,  then  : 

filtering  the  time  section/ 

replacing  the  amplitude  by  one  or  two  bite/ 

filtering  the  new  time  section/ 

If  IENTER  <  3.  then  : 

performing  the  computation  of  the  jump  in  Telocity. 


GE02D. 


This  algorithm  it  directed  by  the  variable  IENTER  which  allows  the 
■aer  entering  the  program  at  different  level  in  the  code  depending  on 
whether  the  naer  is  inputting  data  from  a  data  tape  or  generating  teat  data 
in  line. 

For  IENTER-1  the  execotion  begins  from  the  creation  of  time  traces/ 
that  it  the  scattered  wave  function  will  be  created.  The  wave  function  dne 
to  the  noiae  is  added  to  this  wave  accordingly  to  the  noise-to-signal  ratio 
given  in  the  data.  The  subroutine  NOISE  implements  that  function. 

For  IENTER12  the  wave  function  will  be  filtered  to  provide  high- 
frequency  data  to  the  remaining  code.  This  filtering  part  includes  other 
operations  requested  by  the  formula.  The  modifications  needed  to  generate 
sign-bit  data  are  done  within  the  main  program,  after  the  filtering 
process  (see  the  algorithm  table  II. 2). 

The  inserted  code  generates  a  function  with  sign-bit  amplitude  from 
the  wave  function  created  previously  if  IENTER^l,  or  gathered  on  the  field. 
That  is  the  initial  scattered  wave  is  translated  in  the  Fourier  domain,  to 
be  filtered  with  the  same  filter  as  before.  This  bandlimited  signal  is 
then  studied  in  the  time  domain  and  algn-bit  processing  will  be  applied. 
These  high  frequency  data  will  be  replaced  by  sign-bits  according  to  the 
sign  of  the  amplitudes  of  the  function  if  two  numbers,  -1  and  +1»  are 
need,  or  to  the  value  of  the  amplitude  if  three  numbers,  -1,  0  and  +1,  are 
used  ,  as  follows: 
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ALGORITHM 


II.  MEANS  USED  FOB  HIE  STUD!. 


A.  Ike  coapoter  propria  6E02D. 


The  propria  GE02D  eoapotai  o(x, z)  for  the  2.5  diaemional  aodel.  Thit 
it.  «•  eoniider  the  veve  ipreidiap  in  three  diaemions  bot  ve  coapnte  the 
intepnl  on  the  peophone  irriy  in  tvo  diaemions  under  the  nsnaption  that 
a  is  independent  of  y. 

For  the  pridient  of  a,  called  fi,  ve  pet  the  proceisinp  foranli  in  2.5 
diaemions: 


16nz 


p(i,,)  ~  |  [Be  _  la]  Jdf *{f *F(f )  *exp(-2nif  *2B/c)  * 

|dt’fUJ{,  t)  *exp(2nift) 


(2) 


vhere:  c.  reference  speed/ 

(x. z)  is  •  field  point 

R.  distance  fora  the  field  point  to  the  obaervation  point/ 
f.  frequency/ 

F(f),  filterinp  function/ 
peophone  location/ 

Us($, t),  acatterad  viva/ 

The  chart  in  table  II. 1  ia  an  alporitha  of  the  oripinal  propria 
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* _ 


(1) 


+  «(x,  y,  *)] 

v  c 


where:  -  (x.y.  z)  -  <x.  *)  it  «  field  point/ 

-  v  it  the  welocity  at  the  field  point  (x.y, z)/ 

-  c  it  the  reference  velocity/ 

-  a(x, y, z)  it  the  perturbation  at  the  field  point. 

However  field  data  it  not  foil  bandwidth  data  and  it  cannot  be  simply 
procetted  to  yield  a  eolation  for  a.  Inttead  a  partial  inversion 
consistent  with  the  high  frequency  characterization  of  the  field  data  is 
carried  ont. 

The  Fonrier  transform  of  a  function  contains  information  about  the 
trend  of  a  function  at  the  low  frequencies  or  small  values  of  the  wave 
number  and  information  about  discontinuities  at  high  frequencies  or  large 
values  of  the  wave  number.  Ihus  the  aolution  formula  is  modified  to 
process  only  for  the  discontinuities  in  a(x,y,  z )/  that  is.  the  reflectors 
in  the  subsurface  will  be  located  by  the  peak  amplitudes  of  the  gradient 
function  of  a,  called  0  in  the  following  discussions. 

The  wave  field  received  at  the  geophones  is  the  required  data  to 
compute  the  function  0.  The  latter  behave  like  a  bandlimited  delta 
function  across  each  reflector.  For  true  amplitude  data,  the  peak  amplitude 
data  at  each  reflector  is  in  known  proportion  to  the  reflection  strength. 

The  true  amplitude  data  and  the  sign-bit  data  are  processed  through 
the  Born  inversion  algorithm  to  compute  the  magnitude  of  the  gradient 
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The  purpose  of  this  psper  is  to  show  thst  in  the  sbsence  of  relstive 
trne  amplitude  the  velocity  inversion  algorithm  still  provides  an  image  of 
the  reflecting  snbsorfaces  as  it  was  remarked  in  the  report  of  Bleistein. 
and  Coben  (1984).  Several  noi se-to-signsl  ratios  are  introduced  in  the 
scattered  wave  and  the  noise  effects  are  analysed  in  the  results  of  the 
Born  inversion  algorithm.  Furthermore,  the  reflection  strength  is  compared 
with  the  jump  in  velocity  given  in  the  data  for  each  example  that  will 
have  been  processed  in  the  Born  inversion  algorithm. 

C.  Method  of  approach. 

The  inverse  problem  is  written  for  the  acoustic  wave  equation.  The 
method  was  written  for  a  perturbation  from  a  reference  sound-speed.  A 
nonlinear  integral  equation,  involving  the  product  of  the  unknown 
perturbation  and  the  unknown  wave  field  in  the  earth  was  derived  in  Cohen 
and  Bleistein  (1979)  and  is  also  described  in  Bleistein  (1984).  The 
solution  formula  requires  full  bandwidth  data  which  is  not  available  in  the 
seismic  experiment.  In  fact  the  data  can  be  characterized  as  high 
frequency  data  with  respect  to  most  of  the  length  scales  of  the  seismic 
experiment . 

The  integration  formula  for  the  inversion  uses  Fourier  transforms  to 
yield  the  perturbstion  function  a  defined  by: 


The  papers  which  are  referenced  in  this  literature  review  discussed 
sign-bit  techniques  to  store  the  seismic  data.  Theae  data  were  then 
proceased  through  typical  geophysical  programs.  The  results  were  analyzed 
to  check  how  they  report  information  about  the  aubaurface,  especially  in 
noisy  environments.  This  paper  will  analyze  in  the  same  way  the  reanlts  of 
a  new  mathematical  interpretation  of  the  aeiamic  data.  The  major 
additional  feature  of  this  new  method  is  that  for  true  amplitude  data.  the 
amplitude  of  the  output  is  in  known  relation  to  the  reflection  strength. 
Ihns,  we  can  check  amplitude  accuracy  of  sign-bit  processing  as  well  as 
image  accuracy. 

B. Purpose  and  objectives  of  study. 

I 

We  use  here  a  mathematical  approach  to  seismic  inversion  presented  by 
Cohen  and  Bleistein  (1979)  with  compoter  implementation  in  Bleistein  and 
Cohen  (1984).  The  authors  used  the  Born  approximation  and  asymtotic 
analysis  to  derive  an  algorithm  to  process  seismic  data  and  recover  both 
the  location  of  the  reflectors  in  the  earth  and  to  estimate  the  reflection 
strength  at  each  reflector.  Ibis  approach  bnilds  the  Born  inversion 
algorithm.  An  asymptotic  analysis  of  the  ontpot  of  this  algorithm  was 
made  in  the  same  report.  One  of  the  conclusions  drawn  from  this  analysis 
is  that  errors  in  amplitude  in  the  input  will  degrade  the  estimate  of 
reflection  strength,  however  they  will  not  affect  estimates  of  the  location 


of  the  reflector. 


acceptable/ 

(ii)  high  compression  of  the  data,  aiore  field  processing  controls  are 
possible/ 

(iii)  trne  amplitude  relations  are  maintained/  and  the  procedure  is 
tailormade  for  Vibroseis. 

The  principal  disadvantages  stated  by  the  anthors  are: 

(i)  the  method  is  limited  to  Vibroseis/ 

(ii)  a  heavy  dnty  60Hz  single  phase  power  plant  is  a  possible  source  of 
coherent  noise/ 

(iii)  a  bad  spread  cable  or  array  terminal  can  cause  a  great  loss  of  data/ 

(iv)  the  dynamic  range  is  limited/ 

(v)  a  higher  caliber  field  personnel  is  needed. 

Then  Gimlin  and  Smith  (1977)  and  Smith  and  Gimlin  (1977)  presented  two 
additional  papers  related  to  sign-bit  studies. 

An  analysis  of  sign-bit  recording  for  high  and  low  noise  environments 
was  made  by  Shirley  et  al.  (1983)  for  the  SEG  meeting  in  Las  Vegas,  1983. 
This  approach  to  sign-bit  results  in  a  unified  description  of  the  low 
noise-to-signal  and  high  noise-to-signal  cases.  A  single  expression  that 
is  valid  for  all  random  noise  distributions  shows  the  relationship  between 
input  coherent  energy,  input  random  noise,  and  output  expectation. 

Tuan  Cheng  Lee(1984)  studied  the  amplitude  recovery  with  sign-bit 
data.  Several  examples  including  different  noise-to-signal  ratios  were  run 
and  the  effects  of  the  noise  in  the  amplitude  recovery  were  analyzed. 


conventionally.  Examples  of  actnal  noisy  environnent  conditions  are  shown 
together  with  the  effect  of  superimposing  that  noise  on  simulated  Vibroseis 
data  and  procesaing  it  with  conventional  and  sign-bit  recordings.  They 
simulated  several  types  of  variable  noise  on  the  coapnter  and  noted  that 
they  have  minimal  effect  on  aign-bit  data.  In  the  presence  of  high 
amplitude  background  noise,  such  as  could  be  expected  in  populated  areas, 
on  windy  days,  near  the  surf,  or  in  stormy  seas,  seismograms  resulting  from 
sign-bit  data  are  shown  to  be  significantly  superior  to  those  from 
conventional  floating  point  recording  systems.  Thus,  in  a  typical  noisy 
environment  the  system  recording  sign-bit  only  should  require  significantly 
fewer  memory  modules  and  consequently  reduced  costs  to  acquire  data 
comparable  or  superior  to  those  from  conventional  systems. 

Schoellhorn  et  al .  (197S)  simulated  sign-bit  imalng  in  vibroseis  data 
acquisition  on  the  computer.  A  variety  of  multiple  arrivals  with  differing 
amplitudes  were  evaluated/  these  showed  the  expected  effects  of  amplitude 
ratio  loss,  spectral  loss  when  the  noise-to-signal  ratio  is  less  than  one, 
and  high  impulse  noise  immunity.  However,  field  data  comparisons  of  full 
range  and  sign-bit  summing  showed  only  minor  differences  in  the  final 
stacked  sections. 

An  evaluation  of  sign  bit  seismic  recording  and  its  impact  on 
processing  and  interpretation  was  made  by  Alihilali  et  si.  (197S).  They 
evaluated  the  advantages  and  disadvantages  of  sign-bit  recording.  Some  of 
the  primary  advantages  claimed  by  the  authors  are: 

(i)  smaller,  simpler  spread  cables  and  single  element  geophones  are 


I.  INTRODUCTION 


Sign-bit  data  recording  has  interested  the  seisaic  exploration 
industry  for  aany  years  because  it  provides  an  easy  way  to  greatly  decrease 
the  number  of  channels  in  a  data  acquisition  system.  Sign-bit  digital 
recording  means  that  only  the  sign  of  the  analog  aignal  is  recorded  with 
one  bit.  In  conventional  seisaic  recording  16  to  32  binary  bits  are 
necessary  for  each  sample  point.  So  the  economic  advantages  of  sign-bit 
recording  are  obvious.  It  is  necessary  to  check  that  the  results  obtained 


from  nsual  seismic  processings  are  valid  and  still  provide  good 

„  t 

information.  .  •  ^-1  '  ; ■  jAv* 
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A.  Literature  review. 
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Numerous  published  papers  deal  with  the  data  storage  using  signed  bits 
and  with  the  amplitude  recovery  from  the  information  carried  by  these  data 
on  the  subsurfaces.  Discussions  of  sign-bit  data  recording  in  geophysical 
literature  consist  mainly  of  abstracts  of  papers  at  international  SEG 
meetings.  The  first  paper  was  presented  by  Fort  et  al.  (1973).  Four  papers 
were  presented  in  197S,  namely.  Fort  et  al.  (1975),  Martin  et  al.  (1975), 
Schoellhorn  and  Neale.  (1975),  and  Alihilali  et  al .  (1975)/  Martin  et  al. 
(1975)  claimed  that  the  noise  burst  problems  with  Vibroseis  are  minimised  by 
only  using  the  sign-bit.  Computer  simulation  shows  that  if  only  aign-bit 
data  is  used,  reflection  quality  is  improved  over  that  obtained 


ABSTRACT 


^ Th it  paper  discusses  the  sign-bit  processing  in  the  Born  inversion 


■ethod.  It  shows  on  basic  examples  that  in  the  absence  of  relative  tree 

amplitude,  the  Born  inversion  algorithm  written  by  Bleistein-Cohen  still 
provides  an  image  of  the  snbsnrface.  It  implies  that  all  the  essential 

information  provided  by  s  wave  train  is  contained  in  the  phase  only.  Sign- 

bit  digital  recording  means  that  only  the  sign  of  the  trne  amplitude  signal 
is  recorded  with  one  bit.  In  conventional  seismic  recording.  Id  to  32 
binary  bits  per  sample  point  are  recorded.  The  economic  advantages  of 

sign-bit  acquisition  are  iamediately  obvious.  Complete  amplitude  recovery 
comparable  to  full  gain  recording  can  be  achieved  by  correct  application  of 


sign-bit  techniques. 
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1 )  -  When  the  amplitude  it  positive,  it  it  replaced  by  +1 .  If  it  it 
negative,  its  value  becoaet  -1. 

or 

2) -  The  aaxiaum  aaplitnde  it  stored  in  UMAX.  When  the  staple  value  it 

lets  than  a  certain  percentage  (stored  in  the  variable  PER)  of  ONAX,  it  it 

replaced  by  0.  If  this  value  it  larger, then  the  aaplitnde  replaceaent 

follows  1)  above. 

The  new  signal  or  succession  of  bits  is  then  processed  through  the 

same  filtering  algoritha  nsed  for  the  initial  scattered  wave,  to  eliainate 
low  and  high  frequencies  which  Bight  have  been  introduced  in  the  aign-bit 
processing  . 

To  complete  the  formula  an  integration  loop  ((  in  eq.  2  )  ia  carried 

out  over  a  set  of  processed  traces  for  each  point  in  the  region  of  interest 
in  the  field.  The  data  delimits  this  region  where  it  is  possible  to  image 
the  interface  or  reflector  or  discontinuity  in  the  earth.  For  each  field 
point  (x, z),  the  previous  function  is  divided  by  R*^*  and  multiplied  by  z. 
The  results  of  the  integration  are  scaled  by  a  factor  SCALE,  to  generate  f), 
which  is  then  an  estiaation  of  the  jap  in  velocity. 

The  plot  of  this  function  enables  us  to  study  and  follow  the 
discontinuities  in  the  earth  aodel  which  reflected  the  initial  wave. 
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B.  Description  of  the  date  required  by  the  prograa  GE02D 


The  date  required  to  execnte  the  prograa  GE02D  is  described  in  the 
following  sections. 

The  data  ia  entered  in  free  foraat: 

RECORD  #1 

IENTER :  1  -  SYNTHETIC  DATA  GENERATION. 

2  -  READ  FROM  TAPE7  OF  DECONVOLVED  TRACES. 

3  «=  READ  FROM  TAPE8  OF  PROCESSED  TRACES. 

RECORD  #2 

TMIN :  BEGINNING  OF  TIME  RECORD  MEASURED  FROM  ENTRY  TIME  OF 
SOURCE  PEAK. 

TMAX:  END  OF  TIME  RECORD. 

DT:  SAMPLING  RATE. 

NFFT :  SIZE  OF  FFT. 

RECORD  # 3 

FMINO :  FREQUENCY  AT  WHICH  TAPER  FILTER  BEGINS. 

FMIN1 :  FREQUENCY  AT  WHICH  TAPER  REACHES  ONE  AT  LOW  END. 

FNAX1 :  FREQUENCY  AT  WHICH  TAPER  FILTER  ENDS. 

RECORD  #4 

NGEO:  NUMBER  OF  GEO  PHONES  (TRACES). 

DXSI :  GEOPHONE  SPACING. 

RECORD  #5 

C:  REFERENCE  SPEED. 
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RE COB D  #6 


ZMIN :  SHALLOWEST  Z  (.OE.O)  TO  BE  PROCESSED. 

ZNAX:  DEEPEST  Z  TO  BE  PROCESSED. 

IE:  Z-S PACING  DESIRED. 

RECORD  #7 

XMIN :  LEFTMOST  X  (.GE.O)  TO  BE  PROCESSED. 

XMAX:  RIGHTMOST  X  TO  BE  PROCESSED. 

DX:  X-S PACING  DESIRED. 

RECORD  #8 

NGROUP:  NUMBER  OF  GEOPHONES  TO  BE  HELD  IN  CORE 

RECORD  #9 

MEPRNT :  TRUE  FOR  PRINTOUT  OF  VELOCITX  PROFILE. 

MXTAPE:  TRUE  FOR  MAKING  VELOCITX  PROFILE  TAPE. 

MKPLOT :  TRUE  FOR  MAKING  VELOCITX  PROFILE  PLOT. 

RECORD  #10 

TEETS:  TILT  ANGLE  OF  SXNIHETIC  DATA  PLANES. 

NPLN:  NUMBER  OF  SXNIHETIC  DATA  PLANES. 

NOISIG:  NOISE  TO  SIGNAL  RATIO. 

RECORD  #11 

CS(I) :  SXNIHETIC  DATA  SPEED  ARRAX. 

RECORD  #12 

ZS(I)  :  SXNIHETIC  DATA  DEPTHS  OF  PLANES  AT  X-0 

RECORD  #13 

GOGO:  SET  TRUE  TO  CONTINUE  RUN  VIIH  BAD  DATA. 
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III.  SEVERAL  EXAMPLES. 


Is  every  exeaple  the  synthetic  tiae  section  is  crested  for  s  snbscribed 
snbsorfsce  structure;  it  is  represented  on  the  first  plot  of  esch  following 
section.  Then  these  trsces  sre  processed  through  the  first  pert  of  the 
Born  inversion  slgoritha.  This  process  is  illnstrsted  on  the  plots  nsaed 
'PROCESSED  TINE  SECTION”.  On  the  synthetic  tiae  sections  end  on  the 
processed  tiae  sections,  the  reflectors  sre  locsted  in  tiae/  thst  is  the 
scsttered  functions  peek  st  the  trsvel  tiae  neeesssry  for  the  wsve  to  go 
from  the  geophone  srrsy  to  the  reflectors.  Aftsr  the  trne  saplitnde  dsts 
hsve  been  processed,  the  scsttered  wave  undeletes  only  sronnd  the 
interfsce  location.  On  the  contrary,  after  sign-bit  data  processing  the 
scattered  wave  is  undulating  everywhere  on  the  tiae  section  and  it  is  aore 
difficult  to  locate  in  tiae  the  reflectors.  A  reason  is  that  after  having 
been  filtered  the  synthetic  traces  were  replaced  by  +1  or  -1  all  along  the 
tiae  section  so  thst  the  side  lobes  becoae  as  important  as  the  aain  lobe. 


Then  the  integration  loop  coapletes  the  velocity  coaputation.  At  the 
center  of  the  span  array  of  the  geophones,  thsre  is  just  enoogh  spread  to 
include  all  geophones  that  could  contribute  to  the  answer,  thst  is  to  the 
aaplitude  of  the  velocity.  A  nuaber  of  81  geophones  was  chosen  to  get  a 
good  approxiaation  of  the  velocity.  The  resulting  function  peaks  at  the 
interface  locations,  with  true  nplitude  data  as  well  as  with  sign-bit 


data.  For  each  exsaple,  the  aaplitude  of  the  resulting  function  will  be 


cheeked  at  the  Biddle  of  the  apes  array.  The  recnltt  will  be  diacosaed 
aore  preeiaely  for  each  exaaple. 

A.  The  diacontinoity  ia  a  alngle  horisoatal  plane. 

The  data  aet  for  thia  exaaple  ia  giTen  ia  Appendix  A.  The  plane  ia 
located  2000  ft.  belo*  the  anrface.  The  theoretical  joap  ia  Telocity  doe 
to  thla  diacont inoity  ia  the  earth  ia  900  ft./aec.  The  plota  III.A.l 
illaatratea  theae  data. 

The  reanlta  of  the  filtering  proceaa  are  ah  or  a  on  the  plot  III. A. 2 
when  the  true  aaplitode  data  were  proceaaed  and  on  the  plot  III. A. 4  when 
the  eign-bit  aaplitode  data  were  proeeaaed. 

The  plota  oa  Figorea  III. A. 3  and  III. A. 5  ahow  both  reanlta  of  the 
Telocity  coapntation.  the  foraer  with  troe  aaplitade  and  the  latter  with 
aiga-bit  aaplitode.  Both  calcolatione  locate  the  plane  diacoatiaoity  in  the 
earth  at  the  aaae  depth.  The  aaaerical  reanlta  are  ah  own  ia  table  III. A. 
The  joap  in  Telocity  will  be  read  at  the  abaciaaa  1000ft  froa  the  left  and 
at  the  depth  2000ft.  The  error  in  aaplitode  it  coapoted  froa  the  difference 
between  the  joap  in  Telocity  giTea  in  the  data.  SOOft/aec.,  and  the  joap  in 
Telocity  read  ia  the  reanlta  of  the  Telocity  coapntation. 
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Depth  (feet) 
of  the 
reflector 

2000. 

2000. 

Jaap  (ft/sec) 

at 

1-1 000,  Z -2000 

485. 51 5 

459.31 

error  in 
aaplitnde  (%) 


2.9 


8 


B.  The  discontinuity  is  so  inclined  plane. 

IBs  data  aat  for  this  example  is  given  in  Appendix  B,  where  the  angle 
THE  IS  has  been  replaced  by  10°.  and  ia  illustrated  on  the  plot  III.B.l. 
Before  the  velocity  computation  the  filtering  process  is  applied  on  the 
true  eaplitode  data  giving  the  time  section  illustrated  on  the  Figure 
I11.B.2.  After  having  applied  the  sign-bit  processing  on  the  true 
asiplitude  data  the  same  filtering  process  is  applied  on  the  aign-bit  dsta 
giving  the  results  viewed  on  the  Figure  XII. B. 4. 

The  outputs  of  the  velocity  computation  with  true  amplitude  data 
(Figure  III.B.3)  and  with  sign-bit  data  (Figure  III.B.5)  both  accurately 
locate  the  interface.  The  numerical  results  are  displayed  in  the  table 
III.B/  the  error  is  calculated  from  a  prescribed  jump  in  velocity  which  is 
SOOft/sec. 

C.  Example  with  two  horizontal  interfaces. 

The  data  set  for  this  example  is  given  in  Appendix  C  and  is 
illustrated  on  Figure  III.C.l.  A  first  filter  is  applied  to  the  synthetic 
time  section  and  the  processed  time  section  is  illostrated  on  the  Figure 
III.C.2.  In  a  second  run  the  aign-bit  process  is  applied  to  the  synthetic 
time  section  to  create  aign-bit  data  which  are  then  processed  through  the 
same  filter  and  this  second  processed  time  section  is  illostrated  on  the 
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Figare  III.C.S. 


OUTPUT  OF  THE  VELOCITY  COMPUTATION 
1ITH  ONE  REFLECTOR,  TILTED.  NO  NOISE 


eith 

true  eaplltode 

date 

el  th 

sign-bit  aaplitnde 

data 

Depth  (feet) 
of  the 
reflector 

2000. 

2000. 

mp  (feet/sec) 

at 

;-1300.Z-1770 

464 .795 

440.595 

error  in 
•split ode  (%) 

— di  .64  ' 

-11. » 

H>«  outputs  of  the  Telocity  competition  in  both  rons  tte  illostreted 
on  the  Figure  1II.C.3  end  Figure  II1.C.6.  They  locete  the  interfeces  at  the 
same  depths.  After  both  competitions  the  epper  interface  it  at  the  right 
looation  (1700ft)  and  the  deeper  interface  is  mislocated:  its  depth  is 
2250ft  instead  of  2300ft.  After  having  processed  these  ontpets  in  the 

postprocessing  algorithm  P0STPP,  which  implement!  the  refinements  to  the 
linear  velocity  inversion  theory  (Hagin  and  Cohen,  1984),  both  reflectors 
are  located  approximately  at  the  right  depths  (see  Figures  1II.C.4  and 

III.C.7).  The  numerical  results  are  displayed  in  the  table  1I1.C. 

D.  Example  with  2  tilted  parallel  interfaces. 

The  data  set  for  this  example  is  given  in  Appendix  C,  where  the  angle 
THETS  has  been  replaced  by  10°/  both  interfaces  are  inclined  to  make  10° 
with  the  horixontal.  The  scattered  wave  field  is  on  Figure  III.D.l.  We 
notice  that  the  amplitude  of  the  wave  peaks  shaping  the 
deepest  interfsce  is  smaller  than  the  one  locating  the  upper  subsurface. 

In  the  processed  time  section  (Figure  III. 0.2)  it  is  more  difficult  to 

locate  the  deeper  interface  than  the  upper  one.  The  time  section  with 

sign-bit  amplitude  is  processed  through  the  first  part  of  the  Born 
inversion  algorithm/  the  results  are  represented  on  the  plot  on  Figure 


OUTPUT  OF  THE  VELOCITI  COMPUTATION 
WITH  2  HORIZONTAL  REFLECTORS,  NO  NOISE 


with 

true  aaplitode 

data 

without  |  after 
postprocessing 


with 

sign-bit  aaplitod 

data 

withoot  |  after 
postprocessing 


Depth  (feet) 

11 

of  the: 

1 

l 

-reflector  1 

BBS 

-reflector  2 

mm 

■19 

BS 

1711. 

2320. 


a  (feet/sec) 
-reflector  1 
at  (1000.1700)  498.318 

at  (1000.1711) 
-reflector  2 
at  (1000,2250)  423.738 

at  (1078,  2320) 
at  (1027,2320) 


561 .11 


548.13 


392.393 


407.271 


440.52 


518.80 


error  in 

■  I 

aaplitode  (%) 

1 

-reflector  1 

E mm 

-21.5 

BB 

-reflector  2 

essm 

Dfl 

-18.5 

mm 

TABLE  III.C 


The  outputs  of  the  velocity  conpntation  with  true  aaplitnde  data 
(Figure  I1I.D.3)  and  with  aign-bit  amplitude  data  (Figure  III.D.S)  locate 
the  interfacee  clearly  at  the  ease  depth*.  The  nunerical  reenlt*  are  ahown 
in  the  table  UI.D.  The  error  ia  calculated  after  coaparison  of  the 
coapnted  j asp  with  500ft/aec.  for  the  upper  reflector  and  with  200ft/sec. 
for  the  lower  Interface. 

The  accuracy  of  the  coapnted  j  nap  in  velocity  depends  on  the  distance 
between  reflectors.  Exaaples  were  also  run  with  300ft  and  with  200ft 
between  the  subsurfaces,  showing  that  the  coapnted  jnaps  in  velocity  are 
interrelated:  when  the  distance  gets  saaller,  the  error  in  aaplitnde  gets 
saaller  for  the  lowest  reflector  and  the  error  for  the  nppest  one  gets 
larger. 

Otherwise  the  subsurfaces  were  located  with  the  sane  characteristics  as 
in  the  previous  exanple  III.D/  that  is,  the  velocity  function  peaks  at  the 
right  depth  for  the  upper  subsurface,  and  about  20ft  above  the  prescribed 
depth  for  the  deepest  reflector. 

The  above  reaarks  are  for  true  aaplitnde  data  processing  and  for  sign- 
bit  data  processing. 
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OUTPUT  OF  THE  VELOCITY  COMPUTATION 
WITH  2  REFLECTORS.  TILTED.  NO  NOISE 


vi  th 

troe  aaplitnde 
data 


vi  th 

aiga-bit  aaplitode 
data 


Depth  (feet) 
at  X-Oft 
of  the: 

-reflector  1  1900. 

-reflector  2  2360. 


1900. 

2350. 


Jaap  (ft/aec) 
-reflector  1 
at  (1000,1720) 
-reflector  2 
at  (1000,2180) 


507.7  82 


433.553 


202.642 


425.228 


E.Ihe  scattering  subsurface  it  an  arc  of  circle. 


The  subsurface  it  represented  in  Figure  III.E.l.  The  speed  is  cQ 
above  the  src  of  circle  and  c+  below  it.  The  radios  of  the  circle  is 
aOOOOft  and  its  highest  point  is  located  at  the  depth  hK1000ft  below  the 
center  of  the  array  of  geophones. 

Using  the  data  given  in  Appendix  D,  the  synthetic  tiae  section  (Figure 
III.E.l)  is  created  using  the  following  foraula  which  is  derived  in  the 
paper:  Kirchhoff  field  froa  a  circolar  cylinder  (Bleistein,  1984).  using  the 
TWo-snd-one  half  diaensional  in-plane  wave  propagation  (Bleistein.  1984)/  it 
is  iapleaented  in  the  subroutine  TBACIB  of  the  prograa  GE02D. 


Us(«.£) 


a  B  •exp(2ia(d-a)/c  ) 

_  .  _S _ £_ 

1  b  8n(d-a) 


(3) 


where : 

-  *  is  the  pulsation/ 

-  $*({, 0)  is  the  location  on  the  geophone  array/ 

-  R_  is  the  noraal  reflection  coefficient/ 

-  d  -  ({2  +  (a+h)2)1/2/ 

-  4c(2(d-a))  is  the  geoaetrical  spreading  over  the  two  way  travel  path  froa 
5  to  C  and  back  to  (,  in  the  abscence  of  curvature  effects  of  the 
reflector. 

-  (a/d)^2  is  the  geoaetrical  spreading  doe  to  a  circular  reflector  of 
radina  a  at  depth  h. 

-  2(d-e)/cQ  is  the  two  way  travel  tiae. 

The  processed  tiae  sections  with  true  eaplitode  dsta  and  with  sign-bit 
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data  are  repreaented  on  Figure  III. E . 2  and  on  Figure  III.E. 

The  ootpot  of  the  velocity  coaipntation  with  true 
locatea  the  anbaurface  at  the  right  depth  (Figure  III.E 
computed  with  aign-bit  amplitude  data  (Figure  III.E 
interface  at  the  aaae  depth.  The  numerical  reanlta  are 
table  III.E. 


4  reapectively . 

amplitude  data 
.3).  The  profile 
.5)  locatea  the 
diaplayed  in  the 
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peaks  everywhere  with  large  amplitudes. 


From  this  chapter  IV,  the  following  conclusions  can  be  drawn: 

1.  When  the  noise-to-signal  ratio  is  less  than  or  around  1,  the 

reflectors  are  easily  located  at  the  right  depth  with  sign-bit  data  as  well 

as  with  true  amplitude  data. 

2.  When  the  noise-to-signal  ratio  is  much  larger  than  1,  it  is  more 
difficult  to  locate  the  reflectors  on  the  plot  representing  the  output  of 
the  velocity  computation.  Thia  output  does  not  peak  at  the  right  location 
of  the  interface.  The  error  in  depth  increases  as  the  noise-to-signal 
ratio  gets  larger.  The  same  conclusion  can  be  drawn  with  true  amplitude 
data  and  with  sign-bit  data. 

3.  The  Figure  TV.C.3  illustrates  the  curves  of  the  errors  in  jnmp  in 

velocity  from  true  amplitude  data  processing  and  from  sign-bit  data 
processing  versus  the  noise-to-signal  ratios.  Both  curves  have  the  same 
trend.  When  the  noise-to-signal  ratio  gets  large  the  errors  from  sign-bit 
data  processing  are  smaller  than  those  from  true  amplitude  data  processing/ 
but  when  N/S  is  small  the  former  are  bigger  than  the  latter. 

4.  When  two  bits  are  used  to  store  the  data,  the  reflector  is  easier 

located/  that  is  the  resulting  function  less  ondulates  than  when  sign-bit 
was  osed.  The  accuracy  depends  on  the  threshold. 

5.  Not  filtering  the  sign-bit  data  gives  rise  to  large  errors  in  the 


velocity  computation. 
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OUTPUT  OF  HE  VELOCITS  COMPUTATION 
WITH  ONE  HORIZONTAL  REFLECTOR.  N/S-0.5 
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with 

sign-bit 

data 

with 

threshold 

depth  in  feet 
of  the 
reflector 

2000. 

2000. 

jimp  in  feet/eec. 
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I-lOOOft.ZSOOOft 

368.513 

359.41 

error  in 
amplitude  (%) 

-26.3 

-28.1 

Tabl e  IV. B. 2 
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processing  with  ■  percentsge  eqosl  to  10%  of  the  maximum  signal  amplitude, 
llte  nmeericsl  results  are  displayed  ia  the  table  IV. B. 2. 


The  ease  configuration  is  given  with  a  percentage  PER  equal  to  26%. 
The  synthetic  tiae  section  is  represented  on  Figure  IV.A2.1.  The  processed 
tiae  sections  with  sign-bit  data  and  with  threshold,  on  Figure  TV.A2.3  and 
Figure  IV.B.S  respectively,  are  identically  coaplex.  The  outputs  of  the 
velocity  computation  are  represented  on  Figure  IV.A2.4  when  sign-bit  were 
used  and  on  Figure  IV. B. 6  when  threshold  was  used.  The  amplitudes  of  the 
smiles  are  smaller  in  the  latter  plot  than  in  the  former  one.  The  noise 
influence  has  been  slightly  eliminated  after  threshold  processing  with  a 
percentage  equal  to  26%  of  the  aaxiaoa  signal  aaplitnde.  The  numerical 
results  are  displayed  in  the  table  IV. B. 3. 

C.  The  ontput  resulting  from  noise  added  to  sign-bit  data  is  not 

filtered. 

In  this  section,  we  add  a  noise  with  N/S-0.5  to  the  sign-bit  data. 
Then  the  resulting  function  is  processed  through  the  velocity  computation, 
without  having  previously  been  filtered.  The  plots  IV.C.l  and  IV. C. 2 
lllostrate  the  processed  tiae  section  and  the  output  of  the  velocity 
computation  respectively.  On  the  former  it  is  impossible  to  locate  in  tiae 
the  reflector.  On  the  latter,  the  velocity  fnnction 
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with  with 

•ign-bit  threshold 

date 


depth  in  feet 
of  the 
reflector 


2000. 


2000. 


jump  in  feet/sec. 

•  t 

X«1000f  t.Z -2000ft 


368.513 


370.797 


▼triable  PEB  in  the  data  for  the  prograa  0E02D.  When  the  aaplitudes  are 
below  thia  liait  they  will  be  replaced  by  zero. 

A  eonf iteration  ie  giTtn  in  the  data  in  Appendix  B  with  a  percentage 
PEB  equal  to  5%.  The  aynthetic  tiae  aection  ia  represented  on  Figure 
III.A.l.  In  this  aection.  only  the  coaputatione  using  sign-bit  processing 
are  dicusaed.  The  processed  tiae  sections  with  sign-bit  data  and  with  two 
bits  data.  on  Figore  IV.A2.3  and  Figure  IV.B.l  respectively.  are 
identically  coaplex.  The  outputs  of  the  velocity  computation  are 
represented  on  Figure  IV.A2.4  when  one  bit  was  used  and  on  Figore  IV. B. 2 
when  two  bits  were  used.  The  aaplitudes  of  the  sailes  are  saaller  in  the 
latter  plot  than  in  the  foraer  one. 

The  noise  influence  has  been  slightly  eliainated  after  threshold 
processing  with  a  percentage  equal  to  5%  of  the  aaxiaua  signal  amplitude. 
The  numerical  results  are  displayed  in  the  table  IV.B.l. 

The  same  configuration  is  given  with  a  percentage  PER  equal  to  10%. 
The  synthetic  tiae  section  is  represented  on  Figore  IV.A2.1.  The  processed 
tiae  sections  with  sign-bit  data  and  with  two  bits  data,  on  Figure  IV.A2.3 
and  Figure  IV. B. 3  respectively,  are  identically  coaplex.  The  outputs  of 
the  velocity  computation  are  represented  on  Figure  IV.A2.4  when  two  bits 
ware  used  and  on  Figure  IV. B. 4  when  two  bits  were  used.  The  aaplitudes  of 
the  aailas  are  saaller  in  the  latter  plot  than  in  the  foraer  one. 

The  noise  influence  has  been  slightly  eliainated  after  two  bits 
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OUTPUT  OP  THE  VELOCITY  COMPUTATION 
WITS  (WE  HORIZONTAL  REFLECTOR,  N/S-25.6 


sections  and  the  processed  tiae  sections  sre  very  aessy:  it  is  aore 
difficult  than  before  to  distinguish  the  delta  eaplitnde  corresponding  to 
the  subsorfsce  (see  Figures  IV.A4.1,  IV.A4.2  and  IV.A4.4).  After  the 

velocity  eoapntation,  the  resulting  fnnction  peaks  at  several  depths  with 
about  the  saae  amplitude  at  each  of  thea.  Thu*  it  is  difficult  to  locate 
the  reflector  after  trne  aaplitode  data  processing  as  well  as  after  sign-bit 
data  processing  (see  Figores  IV.A4.3  and  IV.A4.5).  The  nnaerical  results 
are  displayed  in  table  IV. A4. 

the  results  of  executions  with  bigger  noise-to-signal  ratios  shov 
larger  larger  errors  in  the  location  of  the  reflector  as  well  as  in  the 
aaplitnde  of  the  j nap  in  velocity. 

B.  Using  two  bits  to  replace  the  aaplitnde. 

Before  entering  the  sign-bit  processing,  the  created  signal  has 
definite  aaplitndes  everywhere  in  the  tiae  interval/  aost  of  them  are 
saaller  than  the  trigger  value  locating  the  subsurface.  However  all  of 
thea  will  be  replaced  by  *1  and  -1  as  wall  as  the  aaiiaae  veins.  That  is 
why  the  processed  tiae  sections  and  ths  outputs  of  the  velocity  eoapntation 
are  aore  eoaplicated.  That  can  aake  difficult  the  deliaitation  of  the 
interfaces  on  the  plots,  aaong  the  sailss  . 

To  eliainate  this  problea,  ths  aaplitndes  of  the  scattered  wave 
function  will  be  replaced  by  sign-bit  only  when  these  values  are  above  a 
percentage  of  the  aaxian  aaplitode:  this  percentage  is  givsn  by  the 
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OOTPOT  OF  TBE  VELOCITI  COMPUTATION 
WITH  ONE  BOfi E ONTAL  REFLECTOR.  N/S-l. 


with  with 

tree  aaplitode  sign-bit  aaplitode 

data  data 


depth  in  feet 
of  the 
reflector 

2000. 

2000. 

jnap  in  feet/sec. 

at 

X-1000ft.Z«2000ft 

500.235 

359.598 

error  in 
aaplitode  (%) 

-0. 

-28.1 

the  saplitode  of  the  jap  is  saeller 


A2.-N/S-0.S 

Another  exeaple  is  set  with  N/S»0.5.  The  processed  tine  sections  with 
trne  saplitode  dsts  end  with  sign-bit  dsts  ere  represented  on  Fignre 

IV.A2.1  end  on  Figore  IV.A2.3  respectively.  Then  the  velocity  coapotetion 
is  perforaed  with  eech  of  these  processed  tiae  sections/  the  outputs  of 
these  eoapotetions  ere  shown  on  the  Fignre  IV.A2.2  end  on  Fignre  IV.A2.4 
respectively.  The  reflector  is  deliaited  et  the  right  depth  on  eech  of  the 
plots.  The  nnaericel  results  ere  displeyed  in  the  teble  IV. A2. 

A3.-N/S-1.0 

Another  exeaple  is  set  with  N/S-1.0.  The  processed  tiae  sections  with 
true  eaplitnde  dsts  end  with  cign-bit  dsts  ere  represented  on  Fignre 

TV.A3.1  end  on  Fignre  IV. A3. 3  respectively.  Then  the  velocity  coapotetion 
is  perforaed  with  eech  of  these  processed  tiae  sections/  the  ootpnts  of 
these  eoapotetions  ere  shown  on  the  Fignre  TV. A3. 2  end  on  Figore  TV. A3. 4 
respectively.  The  reflector  is  deliaited  et  the  right  depth  on  eech  of  the 
plots.  The  nnaericel  results  ere  displeyed  in  the  teble  IV. A3. 

A4.-N/S-25.6 

Nov  e  bigger  noise  is  crested  with  N/S-25.4.  The  synthetic  tiae 


OUTPUT  OF  TEE  VELOCITI  COMPUTATION 
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wi  th 
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data 
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data 


depth  In  feet 
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reflector 

2000. 
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j map  in  feet/sec. 

at 
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491.697 

384.367 

error  in 
aaplitude  (%) 


-1 .66 


-23.126 


IV.  Effects  of  the  noise  on  the  computation 


A.  Using  one  bit  to  replece  the  saplitnde. 

The  effects  of  the  noise  will  be  studied  in  the  example  of  one 
horizontsl  plane  located  at  2000ft  nnder  the  ground  (see  data  in  Appendix 
B).  When  the  noise-to-signal  ratio  (N/S)  is  not  eqnal  to  zero,  a  wave 
field  is  created  by  the  subroutine  NOISE  and  is  added  to  the  scattered  wave 
field  at  each  sample  point.  The  resulting  wave  field  is  processed  in  the 
same  vay  as  in  the  previous  examples.  Sign-bit  processing  will  be  completed 
with  -*-1  and  -1  as  before. 

A1.-N/S-0.2 

A  first  example  is  set  with  N/S-0.2.  The  processed  time  sections  with 
trne  amplitude  data  and  with  sign-bit  data  are  represented  on  Figure 
IV.A1.1  and  on  Figure  IV.A1.3  respectively.  Then  the  velocity  computation 
is  performed  with  each  of  these  processed  time  sections/  the  outputs  of 
these  computations  are  shown  on  the  Figure  IV.A1.2  and  on  Figure  IV.A1.4 
respectively.  Ihe  reflector  is  delimited  at  the  right  depth  on  each  of  the 
plots.  Ihe  numerical  results  are  displayed  in  the  table  IV. Al. 

The  noise  did  not  affect  the  good  location  of  the  reflector  bnt 
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ODTPDT  OF  IHE  VELOCITI  COMPOTATION 


WITH  AN  ARC  OF 

CIRCLE,  NO  NOISE 

with 

true  amplitude 

date 

with 

aign-bit  amplitude 

data 

depth  in  feet 
it  X-1000 

Z=1000 

1000. 

1000. 

j  amp 

it 

in  feet/sec. 
(1000.  1000) 

4S4 . 

418. 

error  in 
amplitude  in  % 


2 


16.4 


V.  CONCLUSION. 


Based  on  the  conf iguration  examples,  the  following  conclusions  hsve 
been  drawn  concerning  sign-bit  processing  on  ■  linear  backscitter  array: 

(1)  In  the  examples  tested,  the  ontpnt  of  the  velocity  inversion 
algorithm  with  sign-bit  amplitude  processing  is  within  10%  to  20%  of  the 
one  with  trne  amplitude  data  processing. 

(2)  Moreover  the  jamps  in  velocity  computed  in  both  cases  are  of  the 
same  order  as  the  jump  given  in  the  data  for  GE02D,  as  long  as  the  noise-to- 
signsl  ratio  is  not  too  large. 

(3)  In  both  processings  the  noise  does  not  prevent  the  good  location 
of  the  reflecting  subsurface,  as  long  as  the  noi se- to-signal  ratio  is 
rea  sonabl e . 

(4)  For  small  noise-to-signsl  ratios,  the  error  from  sign-bit  data 
processing  is  larger  than  the  error  from  true  amplitude  data  processing. 
Thus  the  noise  bnrst  problems  with  the  Born  inversion  algorithm  are  not 
always  minimized  by  only  using  the  sign-bit,  which  is  in  contrast  with  the 
remark  stated  by  Martin  et  al .  (1975). 
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